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In our present work, five previously proposed sp^ carbon crystals were suggested as silicon allotropes and 
their stabilities, electronic and optical properties were investigated by first-principles method. We find that 
these allotropes with direct or quasi-direct band gaps in range of 1.2-1.6 eV are very suitable for applications 
in thin-film solar cells. They display strong adsorption coefficients in the visible range of the sunlight in 
comparison with diamond silicon. These five silicon allotropes are confirmed possessing positive dynamical 
stability and remarkable themodynamical stability close to that of diamond silicon. Especially, the direct band 
gap M585-silicon possessing energy higher than diamond silicon only 25 meV per atom is expected to be 
experimentally produced for thin-film solar cells. 


PACS numbers: 78.20.Ci, 42.79.Ek, 71.10.-w, 88.40.jj 


Silicon is becoming the fundamental element in our daily 
life. Almost all the semiconductor components in high 
technology of today are silicon-based. As the second most 
aboundant element in earth crust, silicon possesses many 
advantages not only for applications in the semiconductor 
industry. It is also the leading material foundation in the 
important field of photovoltaic energy production. Due to its 
abundance and ability to translate the solar energy, as well 
as the reducing non-renewable traditional energy, silicon is 
believed as the future of human. Theoretically, silicon can 
forms many allotropes due to its ability of sp^ hybridization 
like carbon. But, unfortunately, it is always locked in oxide 
or lie in its diamond like ground state. 

Although the band gap of diamond silicon of about 1.12 
eV— lies in the optimal adsorption range of the sunlight, 
its indirect characteristic seriously affects its adsorption 
efficiency^. Thus, many efforts have been paid on searching 
for viable silicon allotropes with direct or quasi-direct band 
gaps overlaping the sunlight in visible range, for the purpose 
of improving the translating efficiency of the solar energy. For 
example, Botti et al. have proposed some low-energy silicon 
allotropes with quasi-direct band gaps in the range of 0.8-1.5 
eV and possessing strong adsorption coefficients for thin-film 
solar cell applications^. Also, an approach named as Inverse 
Band Structure Design Approach (IBSDA) based on the par¬ 
ticle swarming optimization algorithm (PSO^) was developed 
by Xiang et al. for the purpose of searching for materials with 
expected electronic properties and a quasi-direct band gap 
silicon allotrope Si 2 o-T was predicted^. Based on the IBSDA 
as implemented in PSO, Wang et al. have also found six low 
energy direct and quasi-direct silicon allotropes^. Very recent, 
Lee et al. suggested another approach to search for direct 
band gap silicon allotropes and they successfully discovered 
many direct and quasi-direct band gap allotropes^. Especially, 
a recent experimental progress shows that a cage-like silicon 
allotrope Si 24 with a quasi-direct band gap of about 1.3 eV 
can be synthesized by Na-assistanted two step approach^. 
The crystal structure is confirmed as the same as that of 
the previously suggested CAS-Si^. The authors also show 
that Si 24 possesses strong adsorption coefficient higher than 



FIG. 1: Perspective view of the optimized crystal structues of M585- 
silicon, S-silicon, Z-CACB-silicon, H-silicon and Z-ACA-silicon. 


diamond silicon by first-principles calculation^. All these 
theoretical successful examples make us in confident in using 
computational-alchemy in the engineering of materialsiS. 

In this paper, we predict five low energy silicon allotropes 
(M585-silicon, S-silicon, Z-CACB-silicon, H-silicon and 
Z-ACA-silicon as shown in Fig. 1) based on the crystal 
structures of our previously proposed sp^ carbon phase s*^'^^ 
and find that they are direct or quasi-direct band gap semi¬ 
conductors with strong adsorption coefficients for thin-film 
solar cell application. Such a transplanting idea of crystal 
structures have made a lot of successes in crystal prediction, 
such as zeolite nets to carbon allotropesi^ and carbon crystals 
to silicon allotropes like Bct-C4 silico n^"^:*^ , M4 silico n*"^4^ , 
MIO silicon^ and Cco-SiS^^i^—. Our calculating results show 
that these five silicon allotropes possess positive dynamical 
stability and remarkable themodynamical stability compara¬ 
ble to that of diamond silicon. Especially, the direct band gap 
M585-silicon possessing energy higher than diamond silicon 
only 25 meV per atom is expected to be experimentally 
produced for thin-film solar cells. 

In our present work, we perform the structure opti¬ 
mizations and total energy calculations using the density 
functional theory within local density approximation 
(LDA)22i^ as implemented in Vienna ab initio simulation 
package (VASP)22i22. The interactions between nucleus and 
the 3s^3p^ valence electrons of silicon are described by the 
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TABLE I: Total energies, energy band gaps and fundamental structural information including space group, lattice constants and mass densities 
of diamond silicon, M585-silicon, S-silicon, Z-CACB-silicon, H-silicon and Z-ACA-silicon. 


Items 

Diamond silicon 

M585-silicon 

S-silicon 

Z-CACB-silicon 

H-silicon 

Z-ACA-silicon 

Space group 

FD-3M (No.227) P21/M(No.ll) CMCM (No.63) 

IMMA (No.74) 

PBAM (No.55) PMMN (No.59) 

Lattice a 

5.407 A 

14.64 A 

3.825 A 

3.833 A 

11.832 A 

3.807 A 

Lattice b 

5.407 A 

3.814 A 

17.141 A 

7.289 A 

7.234 A 

7.055 A 

Lattice c 

5.407 A 

6.772 A 

7.399 A 

17.41 A 

3.812 A 

12.002 A 

Mass density 

2.361 Mg/cm^ 

2.311 Mg/cm^ 

2.307 Mg/cm^ 

2.302 Mg/cm^ 

2.284 Mg/cm^ 

2.315 Mg/cm^ 

Relative energy 

0 

25 meV/atom 

42 meV/atom 

70 meV/atom 

61 meV/atom 

78 meV/atom 

Band gap 

1.11 eV 

1.51 eV 

1.41 eV 

1.33 eV 

1.52 eV 

1.29 eV 

Direct band gap 

3.32 eV 

1.51 eV 

1.53 eV 

1.38 eV 

1.63 eV 

1.43 eV 



FIG. 2: Energy per atom as a function of volume per atom for di¬ 
amond silicon, M585-silicon, S-silicon, Z-CACB-silicon, H-silicon 
and Z-ACA-silicon. 


projector augmented wave (PAW) method2i24. A plane-wave 
basis with a cutoff energy of 400 eV is used to expand the 
wave functions and the Brillouin Zone (BZ) sample meshes 
are set to be dense enough (less than 0.21 A“') to ensure 
the accuracy of our calculations. All the crystal structures 
of silicon allotropes considered in our present work are fully 
optimized up to the residual force on every atom less than 
0.005 eV/A. The dynamical stabilities of the hve new silicon 
allotropes were evaluated through simulating their vibrational 
properties by using the PHONON package^^ with the forces 
calculated from VASP. Especially, the hybrid functional 
method (HSE06)S is also considered to accurately calculate 
the electronic and optical properties of the five new silicon 
allotropes. 

As shown in Eig. 1 is the optimized crystal structures 
of M585-silicon, S-silicon, Z-CACB-silicon, H-silicon and 
Z-ACA-silicon. Their relative energy, energy band gaps, as 


well as fundamental structural information including space 
groups, lattice constants, mass densities are summarized in 
Table I , those of diamond silicon are also listed for com¬ 
parison. We can see that the mass density of M585-silicon, 
S-silicon, Z-CACB-silicon, H-silicon and Z-ACA-silicon are 
2.311 Mg/cm^ 2.307 Mg/cm^ 2.302 Mg/cm^ 2.284 Mg/cm^ 
and 2.315 Mg/cm^, respectively. These results indicate that 
M585-silicon, S-silicon, Z-CACB-silicon, H-silicon and 
Z-ACA-silicon are dense close to diamond silicon, which is 
good consistent with their sp^ configurations. 

Erom the total energies as listed in Table I and shown in Eig. 
2, we find that all these hve silicon allotropes are metastable 
with energies higher than that of diamond silicon in 80 meV 
per atom. That is to say, they are themodynamically favorable 
than most of the previously proposed silicon allotropes, such 
as the quasi-direct band gap Si 2 o-T- and Q135ii as well as the 
direct band gap D135i and mC12-Si^. In Eig. 2, we also show 
the total energy per atom as a function of volume per atom for 
these new silicon allotropes (Results about other phase with 
energies higher than diamond silicon up to 80 meV per are 
not shown in it for the considering of simplihcation). Erom 
the quadratic E-V relation near the equilibrium Vo, we can 
see that all of them are themodynamically stable protected by 
obvious energy barriers in their individual local state. Erom 
Table I and Eig. 2, we can easily know that the most stable 
one among these five silicon allotropes is M585-silicon, 
followed by S-silicon and H-silicon. Its energy is just 25 meV 
per atom higher than that of diamond silicon, which indicates 
that it possesses the very high probability to be synthesized in 
experiment. 

To further confirm the dynamical stabilities of these five 
new silicon allotropes, we then investigated their vibrational 
properties through simulation their phonon band structures 
and density of state. Their phonon band structures are shown 
in Eig. 3. According to our results, there is no any imaginary 
frequency appearing in these phonon vibrational spectra and 
no any imaginary models appearing it their phonon density of 
state. Such vibrational properties indicate that these hve new 
silicon allotropes are dynamically stable. 

As mentioned before, direct or quasi-direct band gap char¬ 
acteristics are expected for sunlight adsorption, which can 
avoiding the requirement of thick absorber layer (for indirect 
band gap semiconductors such as diamond silicon) to provide 
necessary phonon moment. Thus, we hrstly care about the 
electronic band structures of these hve promising silicon 
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FIG. 3: Phonon band structures of M585-silicon, S-silicon, Z-CACB-silicon, H-silicon and Z-ACA-silicon. 



FIG. 4: Electronic band structures of M585-silicon, S-silicon, Z-CACB-silicon, H-silicon and Z-ACA-silicon. 


allotropes. Considering that traditional DFT-method always 
underestimate the band gap of semiconductor in comparison 
of the experiment results, we employ the HSE06 method 
in our present work to investigate the electronic properties 
of M585-silicon, S-silicon, Z-CACB-silicon, H-silicon and 
Z-ACA-silicon. The calculated band structures are shown in 
Fig. 4 and corresponding band gaps are summarized in Table 
I. From the band structures we can see that the most stable 
M585-silicon is a typical direct band gap semiconductor with 
a gap of 1.51 eV. S-silicon, Z-CACB-silicon, H-silicon and 
Z-ACA-silicon behave as quasi-direct band gap semiconduc¬ 
tors with band gaps of 1.41 eV, 1.29 eV, 1.33 eV and 1.52 eV, 
respectively. Their corresponding direct band gaps are 1.53 
eV, 1.38 eV, 1.63 eV and 1.43 eV, respectively. These results 
indicate that they are suitable for applications in thin-film 
solar cells as absorbers. 

We then turn our interesting on the absorption abilities of 
these five promising silicon allotropes. Absorption spectra 
of these allotropes are also calculated based on HSF06 
method. The results are shown in Fig. 5 in comparison 
with that of the diamond silicon and the air mass 1.5 solar 
spectral irradiance^^. We can see that all of these five silicon 
allotropes possess strong absorption coefficients better than 
that of diamond silicon in the energy range between 1.5 



Energy (eV) 


FIG. 5: Adsorption spectra of M585-silicon, S-silicon, Z-CACB- 
silicon, H-silicon and Z-ACA-silicon compared to that of diamond 
silicon and the reference air mass 1.5 solar spectral irradiancai^. 


and 3.2 eV. It is well known that indirect band gap diamond 
silicon needs necessary phonon moments to assist its indirect 
band adsorption^. Thus it is always prepared with thickness 
about 100 jum to ensure the adsorption of solar energy below 
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its direct optical gaps22. Our discoveries suggest that these 
five new silicon allotropes are suitable for thin-film solar cell 
applications in views of their direct or quasi-direct band gap 
characteristics and strong adsorption abilities, if they can be 
successfully synthesized in future experiment. 

Usually, silicon is always locked in oxide or lies in its 
ground state of diamond silicon. There is no any graphite-like 
silicon allotrope that can provides abundant pathways to 
many other metastable phases under high pressure condition 
like carbon^. Although diamond silicon can also undergos 
a series of phase transition s^ with the increase of external 
pressure, from diamond to /3-Sn at about 12 Gpa, from pSn 
to orthorhombic Imma phase and then to simple hexagonal at 
about 13-16 Gpa, from simple hexagonal to an orthorhombic 
Cmca phase at about 38 Gpa, from Cmca to a new hexagonal 
close pack at 42 Gpa and finally to face-centered cubic at 78 
Gpa , all these high pressure phases are metallic and non of 
them can keep existing at ambient condition. Fortunately, the 
process of pressure-release of y6-Sn phase provides us many 
opportunities to meet some metastable silicon allotropes^, 
such as R8, BC8 and other unknown phases^ with distorted 
sp^ bondings, depending on the release velocity and temper¬ 
ature condition. We expect further experiments in future can 
find out proper release velocity and temperature conditions 
to synthesize these low energy silicon allotropes based on 
the high pressure y6-Sn phase. On the other hand, other 
experimental methods have also been developed to synthe¬ 
size new metastable silicon allotropes^f^i^. M585-silicon, 
S-silicon, Z-CACB-silicon, H-silicon and Z-ACA-silicon 
are topological similar to the recently synthesized allo-Ge 
phases, which indicate that they may be synthesized by 
similar method^. The recent success of synthesizing of a 
cage-like silicon allotrope Si 24 with a quasi-direct band gap 


of about 1.3 eV by Na-assistanted two step approach^ also 
provides a suitable method for synthesizing the five new 
silicon phases suggested in our present work. Especially, the 
M585-silicon with the lowest energy (about 80 meV lower 
than that of Si 24 ) and cage-like configuration similar to that 
of Si 24 is expected to be synthesized as the first experimental 
direct band gap silicon for thin-film solar cell application. 

In summary, five previously proposed sp^ carbon crystals 
were theoretically suggested as silicon allotropes for the first 
time and their stabilities, electronic and optical properties 
were investigated by first-principles method. We find that 
these allotropes with direct or quasi-direct band gaps in range 
of 1.2-1.6 eV are very suitable for applications in thin-film 
solar cells. They display strong adsorption coefficients in the 
visible range of the sunlight in comparison with diamond 
silicon and remarkable themodynamical stability superior to 
most of the previously proposed direct or quasi-direct band 
gap silicon allotropes such as SUo-T— , Q135i, D135^ and 
mC12-Si^. Especially, the direct band gap M585-silicon 
possessing energy only 25 meV per atom higher than that 
of diamond silicon and about 80 meV lower than that of the 
recently synthesized Si244 is expected to be experimentally 
produced for thin-film solar cells. 
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